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The initial state of a system-bath composite is needed as the input for prediction from any 
quantum evolution equation to describe subsequent system-only reduced dynamics or the noise on 
the system from joint evolution of the system and the bath. The conventional wisdom is to write 
down an uncorrelated state as if the system and the bath were prepared in the absence of each 
other; yet, such a factorized state cannot be the exact description in the presence of system-bath 
interactions. Here, we show how to go beyond the simplistic factorized-state prescription using ideas 
from quantum tomography: We employ the maximum-entropy principle to deduce an initial system- 
bath state consistent with the available information. For the generic case of weak interactions, we 
obtain an explicit formula for the correction to the factorized state. Such a state turns out to have 
little correlation between the system and the bath, which we can quantify using our formula. This 
has implications, in particular, on the subject of subsequent non-completely-positive dynamics of 
the system. Deviation from predictions based on such an almost uncorrelated state is indicative of 
accidental control of hidden degrees of freedom in the bath. 

PACS numbers: 03.65.Wj, 03.67.-a, 03.65.Yz 


I. INTRODUCTION 

Time-evolution problems require the specification of 
initial conditions for the prediction of future behavior. 
In quantum evolution, these initial conditions come in 
the form of the initial state of the quantum system in 
question. For open quantum systems, subsequent dy¬ 
namics involve not just the system of interest, but also 
the bath—the coupling to which gives non-unitary evo¬ 
lution of the system—and the joint initial system-bath 
state is of relevance. How does one write down a rea¬ 
sonable initial system-bath state, given the inability to 
perform tomography on the bath, the defining quality of 
which is that of inaccessibility and uncontrollability? 

Usually, one argues that the system is well isolated 
from the bath: Only when the coupling is weak does this 
split into system and bath make good physical sense. The 
system is then viewed as having been prepared in some 
state ps independently of the bath. The bath, usually 
much larger than the system and hence suffering negli¬ 
gible influence from the system, is taken to be in some 
reference state (e-g-, the thermal state) as if the 

system were absent. This reasoning gives the factorized 
state Ps ® Pb ^he initial system-bath state. One 
arrives at the same state by alternatively imagining that 
the system-bath interaction is zero for time t < 0, and is 
“turned on” only at < = 0. 

Practically, the isolation of the system from the bath 
can never be complete; e.g., one can never decouple from 
the electromagnetic vacuum. One expects the initially 
perfectly uncorrelated situation to be but an approx¬ 
imation, albeit often a good one. As we advance to¬ 
wards more and more sensitive quantum precision mea¬ 
surements, we are less and less justified in ignoring these 
system-bath correlations. The presence of initial corre¬ 


lations also raises interesting questions that have trig¬ 
gered much ongoing discussion regarding the description 
of subsequent system-only dynamics: the lack of a well- 
defined reduced dynamics m-, the connection between 
non-completely-positive (non-CP) dynamics and corre¬ 
lated initial joint system-bath states [IHZ] (see Ref. [4] for 
further references to the many works discussing the con¬ 
sequences of non-CP dynamics in a variety of settings); 
more specifically, the relationship between completely- 
positive (CP) dynamics and initial states with vanishing 
discord [HHin]; attempts to define reduced dynamics in 
the presence of initial system-bath correlations under re¬ 
stricted conditions gmiHIll; etc. The non-CP nature 
of the reduced system-only dynamics resulting from ini¬ 
tial correlations also have close links with non-Markovian 
dynamics of the system, another subject of vigorous re¬ 
cent discussion and creative exploitation of the resulting 
effects (see Refs. [15] and m for comprehensive reviews 
on the subject). An initial correlated system-bath state 
may also have a bearing on conclusions in the discussion 
of topics like quantum fluctuation theorems m- 

In many of these articles, the initial system-bath state 
is unspecified, taken to be some fixed but unknown 
state—unknown because one is unable to perform state 
tomography on the bath. One should, however, be able 
to write down a reasonable initial system-bath state. Af- 
terall, a state is our —the experimenters’—best descrip¬ 
tion of the physical situation. The system-bath state 
should thus express what we know about the system and 
bath, incorporating any prior knowledge of the exper¬ 
imental circumstances, as well as any characterization 
data of the initial system state. Furthermore, it should 
express our lack of knowledge of the precise identity of 
the bath state, stemming from our inability to control the 
microscopic state of the bath, the latter property form- 
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ing the basis for the dichotomy into bath versus system. 
As such, the framework of state estimation [18] provides 
a solid foundation for writing down a joint system-bath 
state. In particular, the approach of maximum-entropy 
(ME) state estimation, well justified in the writings of 
Jaynes [TOHST] and used more recently in quantum state 
tomography , offers a concrete and philosophically 

satisfactory way of representing both our knowledge as 
well as lack thereof for the system-bath composite. 

In this article, we examine this ME route towards a 
reasonable initial system-bath state. The ME approach 
identihes as the best guess for the state, the statistical 
operator that maximizes the entropy in the space con¬ 
strained by our (partial) knowledge from available data, 
with a nonzero entropy arising from the incompleteness 
of the knowledge, i.e., our ignorance. In Sec. II, we set 
the stage for discussion by summarizing the knowledge of 
the system-bath composite obtainable in a typical exper¬ 
imental situation. In Sec. Ill, we explain the use of the 
ME principle for constructing a joint system-bath state. 
We also point out the relationship to the assignment map 
approach to reduced system-only dynamics. In Sec. IV, 
we examine the case of a weak system-bath interaction 
at fixed temperature, the generic experimental situation 
and the case of most relevance in the discussion of ini¬ 
tial system-bath correlations. We obtain a formula for 
the correction to the factorized state, which forms our 
main result. This formula is used to derive a bound 
on the initial system-bath correlations for an arbitrary 
weak interaction Hamiltonian. We discuss implications 
on the subsequent system-only dynamics, and illustrate 
with some examples. In Sec. V, we mention some related 
issues, and conclude in Sec. VI. 


II. OUR KNOWLEDGE 

Controlled quantum experiments begin with the prepa¬ 
ration of the system in some target state at time t = 0, 
before letting it evolve naturally or according to a se¬ 
quence of gates for t > 0. That the system is well- 
initialized into the desired state is verified via tomogra¬ 
phy: The system is repeatedly prepared under the same 
experimental conditions, and the copies are measured us¬ 
ing a chosen tomography strategy; the data are used to 
estimate the prepared state. We phrase this estimate— 
what we know of the system—in terms of estimates of the 
expectation values {oi}^^ of a set of operators 
on the system state. In a typical tomographic setting, 
{Oi} is the set of tomographic outcomes (or POVM) and 
the Oi is the probability that the detector for outcome i 
registers a click. 

The set {Oi} can be complete in that knowing their 
expectation values identihes exactly one density oper¬ 
ator PS for the system state; or it can be incomplete, 
and allows for a restricted set of pss consistent with the 
data. We will assume that {Oi} is complete, although it 
is straightforward to incorporate the case of an incom¬ 


plete set using methods from analyzing incomplete to¬ 
mographic data. We will also assume that we have good 
certainty of our estimates for the o^s. Cases where uncer¬ 
tainty arises because of limited data have been studied in 
the literature and can be incorporated using the concept 
of estimator regions [25ll27j . We will, however, refrain 
from putting in such complications, and only point out 
that one should in general include knowledge about the 
quality of the available information about the system, not 
just what information is available. 

What about our knowledge of the bath? By “bath”, 
we refer to the part of the relevant physical system not 
amenable to microscopic control. We can at best de¬ 
termine only macroscopic properties, like temperature, 
pressure, volume, etc. By “relevant”, we mean that we 
include only those degrees of freedom that couple suffi¬ 
ciently (though weakly) to our controlled system to affect 
its dynamics; the “rest of the world” can be ignored at 
the level of precision detectable in the experiment. We 
possess no microscopic knowledge of the bath, having no 
ability to measure the bath in the same way as we could 
perform tomography on the system. We do, however, as¬ 
sume that we have correctly identified those physically 
relevant degrees of freedom, e.g., the cavity modes that 
are resonant or near-resonant to a few-level atom (the 
idealized system) inside the cavity; the surrounding en¬ 
vironment of nuclear spins around a nitrogen-vacancy- 
centre spin site; etc. 

With this identification also comes the physical model 
for the coupling between the system and the bath, which 
we can phrase in terms of an interaction Hamiltonian 
Hsb- More general coupling, e.g., non-unitary evolution 
via a joint master-equation approach, or a positivity- and 
trace-preserving map, is possible, but for our current dis¬ 
cussion, we will restrict ourselves to the simplest joint 
unitary evolution. Alongside this Hsb, we can also in¬ 
clude a system-only Hamiltonian Hs, and a bath-only 
Hamiltonian Hb- 

In summary, our knowledge of the system-bath com¬ 
posite at time t = 0 can be collected as follows: 

1. information about the microscopic system state ps 
from knowing Oi = tr{psOi}i * = 1; 2,..., AT; 

2. the relevant bath degrees of freedom; 

3. the full Hamiltonian H = Hs + Hb + Hsb', 

4. macroscopic properties of the system and the bath. 

Item I assumes repeatability in the preparation of the 
system state; items 2 and 3 assume a reasonable model 
for the system and bath, based on the experimenter’s un¬ 
derstanding of the physical situation; the last requires the 
laboratory measurement of only macroscopic properties 
of the system-bath composite. 

The above knowledge is the basic starting point for 
any joint system-bath evolution investigation, and hence 
is the right place to begin our construction of a joint 
system-bath state consistent with this information. 
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III. THE MAXIMUM-ENTROPY STATE 

Without the ability to perform tomography on the 
bath to determine its microscopic state, we have only in¬ 
complete information about the overall system-bath state 
at time t = 0. The state that we write down must hence 
reflect this partial ignorance. The use of ME methods in 
such situations follows a long tradition dating back to the 
beginnings of statistical physics. For statistical inference 
problems, Jaynes [HHII] argued lucidly for the viewpoint 
of maximum entropy as the most rational approach in the 
face of incomplete information, and Ref. introduced 
the ME approach to estimate the density operator from 
known expectation values of observables. ME ideas are 
utilized more recently in quantum state estimation with 
informationally incomplete measurements [531121] ■ 

The ME state for the system-bath composite is 

the state that attains 

= =max5(p), (1) 

where S{p) = —trjplogp}, is the von Neumann entropy. 
Here, E is the set of all system-bath states consistent with 
our partial knowledge of the system and the bath. This 
maximization over E can be implemented as an uncon¬ 
strained maximization over all states using the Lagrange- 
multiplier method, so that we instead maximize 

S{p) = S{p) - ^ A*(tr{pOi} - o*) 

i 

-v{tr{p} - 1) - p{tr{G{p)} - g) (2) 

over all system-bath p. Here, AiS are the Lagrange multi¬ 
pliers for the system information, the v term enforces the 
unit-trace constraint for p, and p(tr{G'(p)} — g) is sym¬ 
bolic for the Lagrange multiplier term that accounts for 
any additional knowledge one may have. The ME state 
is obtained from the stationarity condition, i.e., setting 

SS{p) = S{p + 5p) - S{p) (3) 

dp ^logp -k l + i^ + ^A*0, -k I 
to zero for arbitrary infinitesimal 6p. This requires 

C dG \ 

-1 - ^ AiOi - p —j . (4) 

The trace-1 constraint associated with v can be taken 
care of immediately, so that we have the ME state as 


Bayesian corrections to the ME scheme in the determi¬ 
nation of the Lagrange multipliers, when the sample sizes 
are small [5S]. Here, we assume that we have measured 
enough copies so that the A^s can be accurately deter¬ 
mined. 

For each experimental situation, one can view the 
ME strategy as an assignment map that takes 

a system-only state ps to the corresponding system-bath 
state p^^^\ Assignment maps were first discussed in 
Refs. [iHS] as a one-to-one association between system- 
only and system-bath states to permit discussion of re¬ 
duced system-only dynamics. Concretely, with the ME 
assignment map, one can describe the reduced system 
dynamics as a dynamical map 

Ps{t) = ^t{ps), (6) 

with $*(•) = trB{C/(t,0)Al(“^)(-)C/(t,0)^}. 

Here, ps written without any time argument is the initial 
system state verified by tomography on the system, and 
U{t, 0) is the system-bath joint evolution operator for H. 

For any assignment map A, Refs. m considered 
three properties: (a) linear [preserves mixtures, i.e., 
A{pps + {l-p)crs) = pA{ps) + {l-p)A{crs)], (b) positivity 
preserving [Al(ps) > 0 for all ps > 0], and (c) consistent 
[trB{Al(ps)} = Ps]- Pechukas [I] showed that the only 
assignment map that possesses all three properties is the 
factorized map, = PsC)Pb'*^\ for some fixed ref¬ 
erence bath state Pb**^^ . Our ME assignment map , 

generally not of this factorized structure except in spe¬ 
cial circumstances (see below), satisfies (b) and (c), but 
not (a). That it does not preserve mixtures should not 
bother us: When preparing a mixed state ps, one may 
follow a procedure that produces the state directly, or 
one may get ps as a probabilistic mixture of two sepa¬ 
rately prepared states ps,i and ps, 2 - There is no reason 
why the p(“®) states for the two different preparations— 
corresponding to two different physical scenarios—of ps 
should coincide. According to the ME philosophy, the ac¬ 
tual preparation procedure used should form part of the 
prior information entering the constraints. The two dif¬ 
ferent preparations constitute different prior information 
and should hence lead to different pl“®) states. 

As an example of the ME strategy, we consider the 
typical situation where the only controlled macroscopic 
parameter is the temperature T. Such a constraint ap¬ 
pears in Eqs. ([^ and ([^ with p. = /3 = 1/kT {k is the 
Boltzmann constant), G(p) = pH, and g is the average 
energy. The corresponding ME state is 

p(ME) ^ _ PH)/A{. . .}. (7) 



p(“^)=exp(^-A-p^)/tr{...}, (5) 

where A = Here, tr{...}—meaning the trace 

of the numerator when appearing in the denominator 
of a fraction—takes care of the unit-trace constraint as¬ 
sociated with V. Note that there are proposals to use 


Let us examine three cases in this context. 

Case 1. Suppose the bath is in thermal equilibrium, 
and the system is exposed to the bath only at time 
t = 0. This describes many experimental situations 
where the system is prepared in a different part of the 
experimental setup and then moved into the main exper¬ 
imental chamber at time t = 0. The system-bath state— 
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upon introducing the system into the final chamber— 
should hence be the product ps 0 where = 

exp(—/3i7B)/tr{...} is the bath-only thermal state at 
temperature T. Indeed, maximizing entropy subjected 
to the tomography constraints on the system and the 
bath temperature—incorporated using pg{p) —yields the 
expected state. The product structure is clear as S{p) 
involves constraints on system and bath separately, and 
Pg^^ emerges as the usual thermal state from maximiz¬ 
ing entropy for fixed temperature. The appropriate ME 
assignment map for this physical situation is hence of the 
factorized variety, 

A(ps) = PS ® Pb'’^ (8) 


with = p^*^\ 

Case 2. The situation gets more complicated when the 
system is prepared in constant contact with the bath, 
and the combined system-bath is at some temperature 
T. In this case, there is little support for an initially 
uncorrelated system-bath state ps 0 Pb*^^^ , although this 
is commonly used in the literature even for such cases. 
The preparation of the system state can sometimes be 
accomplished very rapidly, but even so, the coupling be¬ 
tween system and bath is always on and can modify the 
outcome of the preparation procedure. Many schemes 
employ slow or long processes that gradually move the 
system population into the desired state, and during this 
transfer, the interaction can be thought of as con¬ 
tinuously keeping the system in thermal equilibrium with 
the bath. In such cases, the coupling to the bath is ac¬ 
counted for through the use of as the initial system- 
bath state for the appropriate T and H. This case is 
examined carefully in the next Section, assuming a weak 
system-bath coupling. 

Case 3. One might also have no macroscopic informa¬ 
tion altogether. In this case, there is no pG{p) term, and 
the ME state is = e“^/tr{...}. Since A is a system- 
only operator, this simplifies to = ps O ;^1 b, 

where dn is the dimension of the bath, ps is the usual 
state one would have deduced from the tomography data, 
with AiS determined only by the constraint summarized 
as trB{p} = Ps- Again, the ME assignment map here is of 
the factorized type with Pb*^^^ = ^1b. The max¬ 

imally mixed state on the bath appears naturally from 
this approach, expressing our complete ignorance of the 
bath state. 


IV. WEAK SYSTEM-BATH COUPLING AT 
FIXED TEMPERATURE 


Let us examine more closely the often-encountered sit¬ 
uation of constant contact between the system and the 
bath, both held at some temperature T (Case 2). We as¬ 
sume a weak coupling between the system and the bath. 
This always-on weak coupling scenario is the case of most 


relevance in the discussion of initial system-bath correla¬ 
tions. Below, we derive an explicit formula for in 

terms of idsBj and explore its consequences. 

Consider a system-bath composite held at fixed tem¬ 
perature T[= l/(/3A:)], evolving with the full Hamiltonian 
H = Hq + Hsb with Hq = Hs + Hb- Hsb is taken 
to be a small perturbation to Hq. The state is 

given by Eq. Q with A determined from the constraint 
trB{p(^^)} = Ps- When Hsb = 0, the ME state is simply 
the factorized state, 

= exp(-Ao - /3ido)/tr{...} = ps O Pb**^ = Alo(ps), 

(9) 

with Ao determined by trBjpg^^^} = exp[—A q — 
/3i7s]/tr{...} = PS so that 

Ao =-logps - (10) 


Here, the natural logarithm log(-) is understood to be 
taken on the support of ps- Aq is a system-only opera¬ 
tor, and is unique up to an additive constant of no con¬ 
sequence: Any constant added to the right-hand-side of 
Eq. (101 is removed by the tr{...} normalization factor 
in Pq ^ . A similar remark applies to A and 

Since Hsb is a weak perturbation to Hq, one expects 
p(ME) close to Po^^^; correspondingly, A is expected 
to be close to Aq. We denote their difference by dA = 
A — Aq, a quantity expected to be small. One can obtain 
an explicit formula for i5A, accurate to linear order in 
/3Hsb ) which we phrase below as a proposition for clarity. 


Proposition. The ME state for weak system-bath cou¬ 
pling specified by the interaction Hamiltonian Hsb o,t 
fixed temperature T = l/{j5k) is 


p(“^) = exp(-Ao - <5A - fiH)M -..}, (11) 

with Aq = — log Ps — fiHs [Eq- ( [To| ) ], and 

SA = -trB{p^“'^/3ifsB}, (12) 


to linear order in j3HsB- 


Proof- Define f{0) = e ‘^/trje for any operator O- 

Let K = Ao+ fiHo- Then, = f{K) and p^^^^ = 

/(A -I- fiH) = f{K 6K), with 6K = dA -I- f3HsB, taken 
as small compared with K. Let — pg^^^ = f{K -\- 

6K) — f{K) = Sf, evaluated as a linear response to a 
variation in K as 


Sf = 


tr{e-^} 


[(5(e ^) -/(A:)tr{(5(e ^)}] (13) 


= -fiK) 


- tT{SKf{K)} 


The second line comes from noting that S{e~^) — 
— dy (<^Ar) from standard perturbation 

theory, and consequently, tr{d(e“^)} = —tr{e“*-diL}. 
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The constraints = ps = that de¬ 

termine A and Aq, respectively, imply trB{(5/} = 0. Note 
that (i) = f{K) can be written as ps ® Pb^^ (ii) 

Pg^^ and K commute; and (iii) A and Aq, and hence dA, 
are system-only operators. Using these facts, straight¬ 
forward simplification of trB{d/} = 0 upon insertion 
of Eq. (13) permits the identification of its solution: 

One requires the relationship dA = —trB{pB*'^/3^^SB}, 


which is Eq. (12). This holds to linear order in /3i?sB, 


due to the linear-response consideration of <5/, which ig¬ 
nores quadratic terms with {SK)^ — {6A -|- /3Hsb)‘^ ^ 

msB?. □ 

Two remarks are in order. First, if the system is pre¬ 
pared in a pure state |^/’s), then, the joint system-bath 
ME state is a tensor-product state by construction. This 
is enforced by the constraint trulp*-^®^} = trB{Po^^^} = 
IV's)('*/'sli via trB{<5/} = 0. This must also be true, accu¬ 
rate to 0{I3Hsb), from direct consid eration of Eq. (11) 
and Eq. (HI- The example in Section [iV C 1| verifies this 
explicitly for the Jaynes-Cummings model. Second, note 
that inserting Eq. (12) for ^A into p(^®) of Eq. (11) yields 
a positive, trace-unity operator for . More precisely, 
one can define an assignment map, 

= exp[logps -6A- P{Hb + i?SB)]/tr{...}, 

(14) 


with SA given by Eq. (12). Then, the ME state for sys¬ 


tem state PS under the weak coupling, fixed temperature 
situation, is simply 


p(ME)=MP(ps) + 0[(/3iPsB)^]. 


(15) 


takes an input system state and outputs 
a system-bath state approximating the ME state. One 
can see explicitly that it is a nonlinear map on ps, ex¬ 
cept in the trivial case when Hsb = 0. For Hsb = 0, 
^0viE)(ps) _ pg (g) pdO _ _ 4 g(pg)^ linearity in 

the input argument is obvious. For the rest of this paper, 
we will use and interchangeably, with the 

understanding that the equivalence is up to 0{/3Hsb)- 


A. Deviation of from uncorrelated state 

That is close to Po'^^^ when Hsb is weak can be 

quantified using the Proposition when the Hamiltonians 
involved are bounded operators: 

Corollary 1. For weak system-bath coupling specified 
by the interaetion Hamiltonian Hsb at fixed temperature 
T=l/{f3k), 

MP(Ps)-pr''^ll<4/3||i?sB||, (16) 


where |jO|| = tr) vO^O > denotes the traee norm. 


Proof. We return to Eq. (13) and rewrite Sf as 
a sum of two terms for the two additive parts 
of SK: 6f = g{6A) -f g{fiHsB), where g{-) = 

-fi^) [/o^ dy eT'^(-) - tr{(-)/(A:)}]. Since the 

trace norm is subadditive and submultiplicative, and 
||/(Ar)|| = 1, we have ||p(0)|| < 2||0|| for any opera¬ 
tor O. This immediately gives ||p(/3iJsB)|| < 2/3|ji?sB||- 
Furthermore, |l5(i5A)|| < 2||dA|| = 2||trB{pB‘'^/3i?SB}|| < 
2^||Hsb||. Thus, ||<5/|| = \\g{SA)+g{l3HsB)\\ < ||5(<5A)|| + 
l!fl(/3^SB)|| < 4/3|ji?sB||, which is Eq. (1^. □ 


This tells us that approximating the initial system- 
bath state as ps C Pb^^ is sufficient as long as we are 
unable to detect effects of order /3||iJsB|i- It is also a 
statement that contains only weak correlations be¬ 

tween system and bath, and hence any effects that rely 
on initial system-bath correlations are small. We see an 
example of this in the next subsection. 


B. Deviation from complete positivity 

An immediate consequence of initial correlations be¬ 
tween system and bath is the possibility of non-CP effects 
in the subsequent system dynamics. A positive map <i> 
is one that takes positive semidefinite operators to pos¬ 
itive semidefinite operators (maps states to states). A 
CP map further requires all extended maps 4) (g) Id to be 
positive, for any d and for Id as the identity map acting 
on an ancillary system of dimension d. The CP property 
ensures that a composite system-ancilla state, generally 
with entanglement between system and ancilla, remains 
positive under 4) that acts on the system only. 

The question of complete positivity comes up natu¬ 
rally when describing system-only dynamics from joint 
system-bath evolution. If the system is initially un¬ 
correlated with the bath, the map that takes 
PS = trB{pSB(t = 0)} to ps{t) = trB{/9SB(t)} = 
trB{U{t, 0)psB(0)U(t, 0)1} is CP. Here, U{t, 0) is the joint 
system-bath evolution operator. When there is initial 
correlation between the system and the bath, as is the 
case when using the ME state as the starting state psB (0), 
may no longer be CP. More precisely, the map will 
be CP if the initial system-bath state has zero discord 
[H]; otherwise, non-CP effects can emerge. Any observa¬ 
tion of non-CP effects can hence be used as a witness for 
the presence of nonzero discord in the initial system-bath 
state [5TI434) . One might also imagine exploiting the ef¬ 
fects that can arise when the CP restriction is lifted, as 
explored in many recent papers. Using our expression for 
^(ME) ^ ggg ggg ggg strong these effects can be. 

A natural “witness” for non-CP dynamics is the trace 
distance between two system states, a measure of distin- 
guishability between the two states. Under a CP map, 
the distinguishability is non-increasing (see, for example, 
[30]'): a violation signals non-CP effects. This fact is used 
to construct a witness for non-Markovian dynamics for 
evolution from time t to an infinitesimally later time t-\-6t 
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[35| : here, we compare only to initial time t = 0, to detect 
the effects of a correlated initial state |36]. The following 
emerges immediately from the Proposition: 

Corollary 2. Consider the situation of weak system- 
hath coupling as in the Proposition. Let be 
a map on system-only states such that ^ti') = 

trB{17(t, 0)^} for some joint evolution U. 

Let ps{t) = ‘ht(ps)- For any two system states ps and 

p's, Zet A(ps,p^;t) = IIpsW - PsWII - IIps - Psll- 

maxA(ps,/Os;i) < 8/3||i?SB|| (17) 

PS.Ps 


for any time t > 0. 


Proof. Consider any two system-only states ps and pg. 


\\psit)-p'sm = \\Mps)-MPs)\\ 

<MP(ps)--4P(p^)|| 

= (ps) - PS o Pb*'^ + Ps O Pb*'^ 

-Ps ® Pb‘’^ + Ps ® Pb**^ - ^r'"^(Ps)ll 
<MP(ps)-pr’'^|| + ||ps-p^|| 
+||(p')r''^-^P(p's)ll 

< 8/3||77sb|| + IIps - Psil- (18) 


Here, (p')o^^^ = -^o(Ps) = Ps C) Pb^\ and we have used 
Eq. (16) in the last line. We thus have that A(ps, pg; t) < 
8/3||i7sB||- This bound is independent of the choice of ps 
and Pg, so Eq. (17) holds. □ 


C. Illustrating examples 

1. The Jaynes-Cummings Hamiltonian 

As a simple example, we work out the ME state for the 
Jaynes-Cummings Hamiltonian, where a spin-1/2 parti¬ 
cle (the system) is coupled to a single mode of light (the 
bath), 


L[ =-huJaCTz+ huja^a +J{a-a^ + a+a). (19) 

2 _^ '—V—" V---/ 

'---" Hb 

Hs 


Hsb 


We denote spin up(down) by tO-)! o-^, i = x,y, z, are the 
Pauli operators, a± = ^{o'x ± io'y) are the spin raising 
((T+ = |t)('H) and lowering (cr_ = | 4 .)(t|) operators, 
and a and are the photon ladder operators. Hsb de¬ 
scribes the absorption and emission of a photon by the 
spin, with the corresponding change of spin state. By an 
appropriate choice of phases in the a and operators, 
or in the spin states, J can be taken as real without any 
loss of generality. Furthermore, for the split into system 
and bath to be sensible, Hsb is a weak correction to the 


free-evolution Hamiltonian, so that J is small compared 
with the energy scales of Hs and Hb. 

The bath-only thermal state is = e“^‘^“^“/tr{...}. 
Since both tr{e“^‘^“^“a} and tr{e“^‘^“^“ai} are zero, 
(5A of Eq. ( [l^ vanishes, so that A = Aq = 
— logPS — \l3huja<yz- The ME state is thus = 

exp [log PS — /3 (Hb -I- HsB)]/tr{...}. A quick check of the 
formula yields = ps®p^^^ = Po'^^^ when Hsb = 0, 

as it should be. 

When Hsb is nonzero, the correction to Pq'^^^ given by 
our ME approach can be worked out explicitly. In this 
case, the ME state can be written as = exp(—AT — 

5K)/tr{- ■ ■}, with K = — logps -\- fiTiuja'^ a and 5K = 
liJ{a-af -\- (T+a). Note that p^^'^ = exp(—Ar)/tr{...} 
in this notation. Using the Bloch-sphere representa¬ 
tion for the state of the two-dimensional system, we can 
write PS = 5(1 -|- s • (t) = e“^''^/tr{...} with A = Ae/, 
where e/ = s/|^, and A = — tanh“^(s). Then, K = 
A • (T -I- phuja^a, where we have dropped an additive con¬ 
stant with no effect on p(^®). Let B be an operator such 
that p(^®) = e^-^pg^^^e^ = e“'®e“'^e®/tr{- • • }; the ex¬ 
istence of such a H is justified below. The operator iden¬ 
tity exp{e~^Ke^) = e~^e^e^ tells us to set e~^Ke^ = 
K -\- 6K = K -\- -\- cr+a). Since /3J is a small 

parameter, this last expression suggests that e® is close 
to the identity, or equivalently, that B is small. Approxi¬ 
mating e^ ~ U-H, we require [K,B] ~ /3 J(CT_a'^-|-cr+a). 
The ansatz B = a+ • aaf -\- a_ • aa -\- u-a-of -\- u+a+a 
satisfies this requirement when 

U± = zLJ/hu! 

and a± = ±i2(J/hw)M±“^(A x /tip), (20) 

where fl± = 5 (l,±i, 0)"^ so that cr± = p± • < 7 , and M± 
are the matrices 


/ Jszphuj —\2\z \2\y \ 

M± = i2Xz ±pnuj -i2X,z ■ (21) 

y —i2Ay i2Aa; Jszjdhoj j 


M± are invertible except in accidental situations 
[det(M) ~ {phuj)^-\-{X terms), which is never zero unless 
A ~ fdhuj]. 

Let us examine two special cases. First, for 
PS = |l, s —and hence A—is zero. Then, B = 
[ J/‘huj){(T-af — cr+a) and p(^®) = pg^^^ — [i3,pg'^^^] = 

p(ME){i _ - 1 ) - a+a(e-^''“ -!)]}> 

to linear order in J/huj. As expected, the ME state 
differs from Pg'^^^ by a term linear in jdJ. If the 
experiment is sensitive to deviations of this size, the 
use of an initially uncorrelated state becomes inap¬ 
propriate. Next, consider a pure state ps = |'0s)('0s|j 
with Bl och vector s of unit length. As mentioned in 
Sec. 


IV 


^^ in Eq. (II) should be a product state 

with trB{p^'^^^} = IV'sTwsI- Solving for Eq. (20), 
taking |s| —>■ I, or equivalently, A —>■ 00 , this gives 
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a± 


p(ME) = , indeed a product state. 

lU 


-prii 

0 

10 “ 

04/3||HsbI| 


0 

0 

0 

2. The central spin model 

10''' 


0 

3 

As a second example, we consider what is known as the 
“central spin model”: a central spin (the system) coupled 

10'^ 

0 

0 



to a number of other spins (the bath). Here, a “spin” can 
be broadly understood as a two-level system or “qubits”, 
e.g., spin-1/2 particles in a bias magnetic field, atoms 

10'^ 

* 

* 

* 

* * * - 


with two energetically relevant levels, etc. The central 
spin model applies in a large variety of physical situa¬ 
tions; see, for example, the review article Ref. [87] . In 
particular, for the development of quantum devices, in 
which the study of the effects of the bath is necessary for 
good control, the central-spin model (and its extensions 
to structured spin baths) has been found to be the ap¬ 
propriate description of noise for many solid-state archi¬ 
tectures; see, for example. Ref. for superconducting 
qubits, Ref. [39] for NV centers, and Ref. [40] for anoma¬ 
lous heating issues in surface traps for ions. These solid- 
state devices are often temperature controlled—e.g., the 
emission/absorption spectral line is temperature sensi¬ 
tive and has to be stabilized—and, through material 
characterization, one often has information on the most 
relevant two-level degrees of freedom forming the spin 
bath. These are the typically available macroscopic in¬ 
formation. 

For simplicity of computation, let us assume all spins 
(one system spin plus N bath spins) have the same prop¬ 
erties, namely, the same energy splitting (of size fno), 
and pairwise interaction (along the quantization axis) 
of the same strength between all spins. The Hamilto¬ 
nian is H = Hs + Hb + Hsb, with Hs = erf, Hb = 

(z) . (z) (j) 1 jj s (2) 

+ 9 ^ and Hsb = 9j2i=i \ 

all measured in units of ^huj, and cri*^ is the Pauli-Z 
operator acting on the ith bath spin. For weak interac¬ 
tion, g <?; 1. One can again apply the same recipe in the 
Proposition to work out the ME state for the fixed tem¬ 
perature situation. Here we focus instead on illustrating 
the two corollaries. 

In Fig. 1, we plot the quantities on the left- and right- 
hand-sides of Eq. (161 in Corollary for fixed g and 
varying the number of bath spins N. One sees that as 
N increases, the difference between ME state and 
increases. Meanwhile, the bound given by Eq. (161 in¬ 
creases faster. This is not surprising, since our bound 
is a uniform bound that holds for all situations. Figure 
2 shows the same quantities as in Fig. 1, but now with 
fixed iV = 1 and varying g. As the g decreases, we see 
the two states approaching each other. Note that the 
blue circles and red asterisks differ by a g-independent 
constant order of magnitude, consistent with a scaling of 
both quantities by the same power of ||/3i7sB||- 

As an illustration of Corollary 2, we study a similar ex¬ 
ample of a central spin coupled to a bath spin, but with a 


10 “ 


4 

N 


FIG. 1. A central spin coupled to N bath spins: The left-(red 
asterisks) and right-hand-sides (blue circles) of Eq. (16) as a 
function of the number N of bath spins in the central spin 
model, g is fixed at 10~^. ps is a random system state; the 
features of the graph do not vary much with different choices 
of Ps. 


10^ 


10'“ 


10 


10"' 


10'^ 


*\\Ar\ps)-pr\\ 

0 4/3||Hsb|I 

( 

0 

\ 

0 

a 

° * 
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10 ‘‘ 


10"' 


10"' 


10 '- 


10 


10 '- 


FIG. 2. A central spin coupled to one bath spin: The left- 
(red asterisks) and right-(blue circles) hand side of Eq. (16) 
as a function of the interaction strength g in the central spin 
model. As in Eig.j^ ps is a random system state; the features 
of the graph do not vary much with different choices of ps- 


randomly chosen Hamiltonian H. Fig. shows the typi¬ 
cal behavior of A(ps, Ps ', 0 ^ randomly chosen pair of 

initial system states. Eq. 0 gives the maximal viola¬ 
tion of the non-increasing property of distinguishability 
under a CP map for any two initial system states. It is 
a precise statement of how initial system-bath correla¬ 
tions in the weak coupling situation have only weak con¬ 
sequences: One needs to detect effects of order /I|jiJsB|| 
to see a deviation from dynamics starting from a factor¬ 
ized initial state. 


V. DISCUSSION 

The system tomography process reliably verifies that 
the system is prepared in a certain state ps- For each ps. 




















X 10"^ 



time t 


FIG. 3. A central spin coupled to one bath spin: The full 
Hamiltonian H is randomly chosen with the requirement that 
Hsb is 1% in norm compared to Hs and Hb; the two initial 
system states are also randomly chosen. The plot shows a 
typical graph for A{ps, Ps;t). One sees that A hardly rises 
above 0 (t in units of the maximum eigenenergy of the system 
divided by h). While Eq. (171 provides an upper bound for 
A, we note that this bound is not anywhere close to being 
achieved for any of the cases in our numerical simulation of 
3000 randomly chosen Hamiltonians, each with 100 pairs of 
system initial states. 


there are infinitely many system-bath states p consistent 
with trel/o} = ps- Macroscopic constraints on the sys¬ 
tem and the bath further restrict p, giving S as the set 
of system-bath states consistent with the available infor¬ 
mation. The inability to measure microscopic features of 
the bath translates into our inability to single out—based 
only on experimental data—a particular p from S. In the 
preceding sections, we employed the additional principle 
of maximum entropy to identify a unique p from S, as a 
quantitative statement of our lack of knowledge. There, 
any knowledge of the system-bath state is accounted for 
via the use of Lagrange multipliers. 

We note that our knowledge and partial ignorance of 
the system-bath state can alternatively be incorporated 
using Bayesian methods, provided one knows how to 
choose the underlying weight properly. Concretely, one 
reports the Bayesian mean well known from Bayesian 
statistical inference: the average over all p in E, p = 
/j,(dp)p. Here, (dp) is a normalized weight density— 
the prior—over E, chosen to represent our available in¬ 
formation. For example, consider Case 3 of Sec. |III| where 
no information other than the system tomography data 
is available, so that E comprises all states such that 
ti'nlp} = Ps- A plausible expression of this complete 
ignorance of the bath state is the Haar average construc¬ 
tion, p = /(df7B)(ls o t/B)|V')sB(V’l(ls o ul), where 
f (dC/e) averages over unitaries on the bath according to 
the Haar measure, with |'0)sb a fiducial state such that 
trB{|l/’)sB(V'|} = Ps- This gives p = ps^Ib/^b, identical 
to that obtained earlier using the ME reasoning. How¬ 
ever, we note that a different choice of “uniform” weight 


over the bath states generally yields a p different from 
the corresponding p(^®), and care should be taken in the 
choice of (dp) for an apt expression of our ignorance of the 
bath state. We note also the work in Refs. gniig based 
on minimizing the relative entropy to a state r which 
is treated as the prior; their procedure yields the ME 
state for the typical experimental situation of system- 
bath composite held at fixed temperature, provided r is 
the joint system-bath thermal state. 

That the initial correlations between the system and 
bath are weak, as in Sec. IV A is a consequence of the 
weak Hsb- It is tempting to extend this to the case of 
strong system-bath coupling. Certainly, there are four 
conditions under which the correlations between system 
and bath will clearly remain weak even for strong Hsb- 

(a) the state ps is nearly pure, so that any correlations 
with the bath state must necessarily be weak, regardless 
of Hsb', (b) the system is prepared separately from the 
bath, so that the product ps® Pb**^ is the appropriate 
initial state (Case 1 of Sec. HI); (c) there is no knowledge 
other than the system tomographic constraint, such that 
the ME state is the product state ps 0 ^1 b (Case 3 of 
Sec. HI); (d) there is a temperature constraint on the 
system and bath, but the temperature is high enough 
such that many system-bath states are involved in the 
thermal mixture that follows from the ME state. 

When none of the four conditions above hold, and Hsb 
is strong, the weak system-bath correlations conclusion is 
unlikely to be the right one. In this case, in the absence 
of the system tomographic constraint as enforced by A, 
the system-bath state will have a large weight on the 
ground state of the full Hamiltonian H, which includes 
the effects of the strong Hsb- With the state prepa¬ 
ration constraint [but not satisfying Conditions (a) or 

(b) above], the system-bath state can still generally have 
strong correlations. One way of treating this would be 
that of Ref. [2S] and others, of treating the state prepa¬ 
ration as an operation done on the possibly strongly cor¬ 
related joint thermal state. However, in such a situation, 
the split into system and bath is called into question. For 
quantum information processing tasks, in particular, one 
cannot hope for good control if one chooses to work with 
a system strongly coupled to an uncontrollable part. 

So far, we have been considering a single joint system- 
bath state consistent with our incomplete knowledge. 
Such “point estimators” are useful whenever one prefers 
a single “best guess”, for use in future calculations. One 
might instead be interested in properties of E itself. For 
example, suppose we only have knowledge about the 
system state ps, and one is concerned with the ques¬ 
tion: How far can any state in the set E be from a 
tensor-product state? In this case, we can make use of 
the conditional entropy 5(B|S) = 5(B,S) — 5(S), where 
5(B,S) = 5 (psb) is the joint system-bath entropy. The 
conditional entropy captures how much of the bath state 
we know, given ps - A tensor-product state ps 0 Pb gives 
iS(B|S) = S{pb), which is bounded by the maximal value 
of log^B (attained by the state ps 0 1b/^b G E), and 
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the minimal value of 0 (attained by ps ® Pb € S when pB 
is pure). Any negative value of 5(B|S) hence indicates 
correlations (classical or quantum) between the system 
and bath. The lower bound is attained by the most cor¬ 
related state in S, namely a purification of ps, so that 
minps3gE5(B|S) = -5(ps). 

In fact, the maximum entropy state that we look for 
in Sec. |III| also maximizes the conditional entropy 5(B|S) 
as the entropy of the system state is fixed by the tomog¬ 
raphy. It thus comes as little surprise that the maximum 
entropy state is uncorrelated or almost uncorrelated be¬ 
tween the system and the bath. On the other hand, if 
one has at hand a situation requiring a conservative es¬ 
timate of the maximal amount of correlation consistent 
with the given data—e.g., in the case of security analysis 
of quantum key distribution protocols—then a purifica¬ 
tion |'0)sB of Ps in E can be used. However, one should 
remember that such a pure state can only serve to derive 
bounds for the “worst-case scenario” and should not be 
misconstrued as a description of the actual system-bath 
state. The lack of control on the bath should translate 
into the lack of our ability to single out a particular |^/')sb- 

VI. CONCLUSION 

The initial joint system-bath state is not a myste¬ 
rious quantity that is unknown or arbitrary. Instead, 
one should take proper account of our knowledge of the 
system-bath composite, encompassing all measurements 
already taken to characterize the experimental circum¬ 


stances. This knowledge gives us a good handle on the 
system-bath state through the use of ME ideas from 
quantum tomography. In this article, we derived an ex¬ 
pression for the ME state in the case of weak system- 
bath coupling, and showed that it is close to a factorized 
state. The rationale for the ME strategy lies with the by¬ 
definition uncontrollable nature of the bath: Maximizing 
the entropy is a quantitative statement of that uncontrol¬ 
lability. Any observed deviation from the predictions of 
the ME joint state should thus be viewed as indicative of 
the presence of some microscopically controlled, repeat¬ 
edly identically prepared degrees of freedom in the bath 
that one failed to identify. Discovery of these hidden con¬ 
trolled degrees of freedom would allow for the possibility 
of exploiting them for better and more creative system 
control. Otherwise, the ME state provides a good de¬ 
scription of the initial system-bath situation, with weak 
correlations in typical situations of experimental interest. 
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